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Using high temperature series we calculate temperature derivatives of the spin-spin and density- 
density correlation functions to investigate the low energy spin and charge excitations of the two- 
dimensional t-J model. We find that the temperature derivatives indicate different momentum 
dependences for the low energy spin and charge excitations. By comparing short distance density- 
density correlation functions with those of spinless fermions and hard core bosons, we find that the 
t-J model results are intermediate between the two cases, being closer to those of hard core bosons. 
The implications of these results for superconductivity are discussed. 



The nature of the ground state and low energy exci- 
tations for two-dimensional strongly correlated electrons 
doped slightly away from half-filling is of considerable 
interest for understanding the properties of high temper- 
ature superconductors Q . The t-J model on a square lat- 
tice is a widely studied model used to investigate these 
problems. While the properties of a single hole intro- 
duced into an antiferromagnet are well understood 
how to extend these results to a finite density of holes 
remains a subject of much current research S. 

For conventional metals, the electron spectral function 
is the simplest way to investigate the energy and mo- 
mentum dependence of the single particle excitations [El . 
With high temperature series we cannot easily calculate 
the spectral function of the 2D t-J model directly. From 
the high temperature series for the momentum distri- 
bution 7ik we calculated || the temperature derivative 
<ink / dT which we used as a proxy for the momentum de- 
pendence of the low energy part of the spectral function. 
Our results for dn^/dT showed that the low energy ex- 
citations of the 2D t-J model are spread throughout the 
Brillouin zone and are in general not conventional quasi- 
particles. A consequence of this result is that dn^/dT 
does not completely determine the momentum depen- 
dence of the low energy elementary excitations in the 2D 
t-J model. To further investigate the nature of the low 
energy elementary excitations, we have extended our cal- 
culations to the equal time spin-spin and density-density 
correlation functions, 5*(q) and N(q) respectively, and 
their temperature derivatives dS(q)/dT and dN(q)/dT. 

We have calculated high temperature series for 5(q) 
and iV(q) of the 2D t-J model to twelfth order in inverse 
temperature f3 = 1/fcsT. Our calculations extend previ- 
ous series calculations ||[7| for the correlation functions 
of the t-J model. The t-J Hamiltonian is given by 
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where the sums are over pairs of nearest neighbor 
sites and the projection operators P eliminate from the 
Hilbert space states with doubly occupied sites. The def- 
initions of the spin-spin and density-density correlation 
functions are 

S(q)=5>^ r (S ^} (2) 

r 

N( q )=J2^ r (An An r ), (3) 



where S* = \ £ Q/3 4 a a z af3 c r p and An r = £ CT c\ a c Ta - 
n. The series are extrapolated to T = 0.2J by Pade 
approximants and a ratio technique used previously || 
for 7ik- 

To interpret our results for the t-J model we exam- 
ine dN(c\)/dT for the tight-binding and spinless fermion 
models. For these non- interacting models dN(q)/dT is 
given by 
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where g = 2 for tight binding and g = 1 for spin- 
less fermions. The properties of dN(q)/dT for the non- 
interacting models are determined by the convolution of 
nk with drik+q/ dT . At low temperatures, due to Fermi 
statistics nk is large only inside the Fermi surface and 
dnk+q/ dT is negative just inside the Fermi surface, pos- 
itive just outside and close to zero elsewhere. The con- 
volution will then give a significant contribution to the 
integral in Eq. 4 when q is such that only one of the pos- 
itive or negative parts of (ink+ q / dT overlaps n^. This 
occurs for q « and q « 2kp, as demonstrated in Fig. 1 
where we plot dN(q)/dT for the tight binding and spin- 
less fermion models. The main features of these plots are 
a large positive spike at q « and a smaller but more 
extended negative dip located at q w 2kp. The shape 
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of the 2kp line depends on the nature of the Fermi sur- 
face (hole like or electron like) but in both cases 2kp is 
a continuous curve in the Brillouin zone. 



at T = 0.2 J and T = 0.4 J we calculate A7V(q)/AT 
and AS f (q)/AT as approximations for the temperature 
derivatives at T — 0.3J. 




FIG. 1. (color) Full Brillouin zone plots of dN(q)/dT for 
left, the tight binding model and right, the spinless fermion 
model. The parameters for both models are n = 0.8 and 
T = 0.08t. Large, positive central spikes show up as red and 
dark blue, while the negative features at 2kp are orange and 
yellow. 

For the t-J series calculations the q«0 (long range) 
parts of the correlation functions have the least accu- 
racy, while we expect the correlation functions at larger 
wavevectors (short range) to be well determined. Conse- 
quently, we concentrate in our analysis on the locations 
and properties of the "2kp" features in the t-J model 
correlation functions. Using the non- interacting models 
as guides, we search for the "2]s.f" features in the t-J 
correlation functions by looking for the largest negative 
values of dJV(q) /dT and dS(q) fdT. 




FIG. 2. (color) Three dimensional full Brillouin zone plots 
of correlation functions for the t-J model. Top, S(q) and 
bottom, iV(q). The parameters for both plots are n = 0.8, 
J/t = 0.4 and T = 0.2J. 

Results for the t-J model iV(q) and S(q) with elec- 
tron density n — 0.8, J/t = 0.4 and T = 0.2 J are shown 
in Fig. 2. Our data are in good agreement with previ- 
ous calculations at higher temperatures (^J^]. The Bril- 
louin zone sums of the correlation functions agree with 
their respective sum rules to within 0.5%. Using data 




FIG. 3. (color) Two views of the full Brillouin zone plot 
of AS(q)/AT in units of J" 1 for the t-J model. The data 
and color coding are the same for both plots, with parameters 
n = 0.8, J/t = 0.4 and T = 0.3J. 
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FIG. 4. (color) Two views of the full Brillouin zone plot of 
AiV(q)/AT in units of J -1 for the t-J model. The central 
peak with a maximum value of 0.08J - has been truncated 
at 0.01J -1 in the lower plot to better show the behavior at 
large wavevectors. The data and color coding are the same for 
both plots, with parameters n = 0.8, J/t = 0.4 and T = 0.3 J. 

Results for AS(q)/AT are shown in Fig. 3. The 
peak at q « is considerably broader than for the non- 
interacting models, in agreement with the temperature 
dependence observed in previous calculations ||. The 
only part of the Brillouin zone where A5(q)/ AT is neg- 
ative is a roughly circular region of approximate radius 
0.67T centered on (tt, it). We note that even at the lowest 
temperature accesible to us, the spin correlations are still 
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peaked at (tt,tt). The correlation length, around (%, tt) 
appears to be saturating or perhaps even decreasing with 
lowering T This could be taken as an indication 
for incommensurate correlations at still lower tempera- 
tures, as the peak at (tt, tt) splits into several distinct 
peaks. Experimentally, for the high temperature super- 
conducting materials, incommensurate spin-correlations 
only arise below about 100 K 

The negative feature in AS'(q)/AT does not form 
a closed curve in the Brillouin zone. In particular, 
A5(q)/AT remains positive from (0, 0) to (tt, 0), with 
no indication of a Fermi surface in the low energy spin 
excitations along this line. Interpreting part of the neg- 
ative feature in A5*(q)/AT as due to an underlying mo- 
mentum distribution of itinerant spin degrees of freedom 
gives disconnected arcs of low energy spin excitations cen- 
tered near (7r/2, 7r/2) and extending perpendicular to 
the zone diagonals. These features are located near the 
peaks observed J5| in dn^/dT and |Vk^k|, consistent with 
the strongest features in being due to an underlying 
spinon Fermi surface. 
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FIG. 5. Top, N(n, 0) and bottom, N(ir, tt) versus tem- 
perature. Solid line: spinless fermions, squares: hard core 
bosons and triangles: t-J model. The density is n = 0.8 for 
all of the data sets and J/t = 0.4 for the t-J model. 

Results for AA(q)/AT are shown in Fig. 4. The peak 
for q « is much sharper than for AS'(q)/AT and more 
like the non-interacting models. The negative feature in 
AA(q)/AT does make a closed curve in the Brillouin 



zone, in contrast to AS(q)/AT. The location and shape 
of the negative feature in AA(q)/AT are similar to the 
2kp line in dN(q)/dT of the spinless fermion model at 
the same density. The similarity extends to having the 
strongest negative feature in both models near (tt, 0). In 
the spinless fermion model this is due to parts of the 2kp 
curve overlapping after being translated back into the 
first zone. The momentum width of the negative feature 
in AiV(q)/AT for the t-J model is considerably broader 
than the 2kp line for spinless fermions and this width 
is not temperature dependent down to T — 0.2 J. Inter- 
preting the negative feature in AN(q) / AT for the t-J 
model as arising from an underlying momentum distri- 
bution for the charge degrees of freedom gives low energy 
charge excitations smeared out over a range of momenta 
near k^? of spinless fermions at the same density. For 
n = 0.8 the charge excitations are centered on the zone 
diagonals and away from (tt, 0). 

The low temperature momentum dependence of N(q) 
for the t-J model, while in general similar to spin- 
less fermions fij, has important differences from spinless 
fermions at large momenta. As shown in Fig. 2, N(q) for 
the t-J model near (tt, tt) and (tt, 0) is slightly smaller 
than 0.2, the value of N(q) for spinless fermions when 
q > 2kp. This means that holes in the t-J model have a 
greater tendency to be nearest or next-nearest neighbors 
than in the spinless fermion model, with the effect largest 
for nearest neighbors. This tendency has also been ob- 
served in exact diagonalization Q and Green's function 
Monte Carlo calculations , in both of which the effect 
is much larger than observed in the series calculation. 
This enhancement is probably due to finite size effects, 
though more work is needed to fully understand the dif- 
ferences between the series calculations and the Green's 
function Monte Carlo and exact diagonalization calcula- 
tions. 
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FIG. 6. Real space, nearest neighbor hole-hole correlation 
functions. Solid line: spinless fermions, squares: hard core 
bosons and triangles: t-J model. The density is n = 0.8 for 
all three data sets and J/t — 0.4 for the t-J model. 

Fig. 5 shows the temperature dependence of iV(q) at 
(tt, 0) and (tt, tt) for the t-J, spinless fermion and hard 
core boson models. All of these models have an infinite 
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on-site repulsion which sets the overall scale for N(q). 
The hard core boson results shown in Figs. 5 and 6 are 
derived from a twelfth order high temperature series for 
AT(q) of hard core bosons. The ground states of spin- 
less fermions and hard core bosons are known: Fermi sea 
with a well defined Fermi surface for spinless fermions 
and a superfluid for hard core bosons. The data for the 
t-J model lie between these two cases, leaving open the 
possibility that the t-J model has a superfluid ground 
state. The real space nearest neighbor hole-hole corre- 
lation function (#0^(1,0)) further supports this behavior. 
Fig. 6 shows the temperature dependence of (<$o^(i,o)) 
for the t-J , spinless fermion and hard core boson models. 
Again, at low temperatures the t-J data is between the 
spinless fermion and hard core boson results. 

The temperature dependences shown in Figs. 5 and 
6 are also interesting. At high temperatures all three 
models have similar values, with the t-J data very close 
to spinless fermions. As the temperature is lowered be- 
low T ^ 1.5J (for J/t = 0.4) the t-J data deviates from 
spinless fermions towards hard core bosons. This tem- 
perature scale is too high to be due to coherent spin fluc- 
tuations. Also, the t-J results are only weakly dependent 
on J/t for J/t SS 0.5 and persist to J/t = 0. This shows 
that the enhancement of (<W(i,o)) m the t-J model rel- 
ative to spinless fermions is due to the presence of two 
spin species, but not due to a direct spin interaction. 

For larger J/t, in the phase separation parameter 
regime, the nearest neighbor density correlation grows 
rapidly as the temperature is lowered. This latter in- 
crease is clearly due to an attraction between the holes 
mediated by J. In contrast to this, the high tempera- 
ture devitaion from spinless fermions towards hard core 
bosons appears to be a statistical effect. This is, per- 
haps, analogous to the Hanbury-Brown and Twiss cor- 
relation well known for quantum particles [[Il| , [I2| . Thus 
it appears that strong on-site repulsion tends to make 
the charge degrees of freedom in the t-J model similar to 
hard core bosons at a high temperature scale JT3|. This 
would suggest a superfluid ground state for the model, 
where the spin degrees of freedom merely help to choose 
the symmetry of the superfluid state. For antiferromag- 
netic spin correlations d-wave might be favored, while for 
ferromagnetic spin correlations p-wave symmetry could 
be favored. These ideas find some support in two-hole 
calculations (l4| , but require considerable further inves- 
tigation. 

In conclusion, from scries calculations we find that 
AS(q)/AT and AA(q)/AT for the 2D t-J model 
are quite different, giving further support to non- 
quasiparticle elementary excitations in 2D strongly corre- 
lated systems. Interpreting these results as due to itiner- 
ant degrees of freedom with underlying momentum dis- 
tributions gives low energy spin and charge excitations 
near the zone diagonals, but with different momentum 
dependences, and absent near (it, 0). The charge correla- 
tions show considerable similarity with hard-core bosons, 
which is suggestive of a superfluid ground state. 
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